Generic features of the phase transition in cold and dense quark matter 
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We investigate the phase transition in cold and dense quark matter in an intuitive way that shares 
common features of the effective model approaches. We first express the quasi-particle contribution 
to the thermodynamic potential with the dynamical mass M and then discuss how we can understand 
the possible first-order phase transition with and without the vector interaction from the saturation 
curve on the plane of the energy per particle and the density. We next extend our analysis including 
inhomogeneity and discuss the relation between the order of the phase transition and the saturation 
curve. We emphasize that the saturation curve is useful to infer qualitative nature of the phase 
transition even without knowing the explicit solution of the ground state. 

PACS numbers: 21.65.Qr, 12.38.Mh, 25.75.Nq 



I. INTRODUCTION 

The quest for the phase diagram of strongly interact- 
ing matter out of quarks and gluons (i.e. matter de- 
scribed by Quantum Chromodynamics - QCD) is one 
of the most challenging problems in modern theoreti- 
cal and experimental physics. There are many specula- 
tions on the QCD phase diagram from theory such as the 
color-superconducting phase [1], the quarkyonic state [2] 
and the triple-point-like structure [3], the QCD critical 
point [4], and so on (see Refs. [5, 6] for comprehensive 
reviews). Available experimental information [7] is, how- 
ever, too limited to constrain uncertainties on those spec- 
ulative possibilities (see Ref. [8] for an attempt and also 
Refs. [9, 10] for physical interpretations). Among others, 
the QCD critical point search is vigorously ongoing in the 
present and future experimental facilities as well as in the 
first-principle calculation of the lattice-QCD simulation. 

The QCD critical point would be, if discovered, a land- 
mark for our understanding on QCD matter. In an 
infinite-volume system at equilibrium, the fluctuations 
are expected to show critical behavior and thus the criti- 
cality would serve as experimental signatures [11]. There 
are many theoretical proposals and experimental data 
taken from the beam-energy scan program at Relativis- 
tic Heavy Ion Collider (RHIC) in the Brookhaven Na- 
tional Laboratory. It is an urgent question to make clear 
whether the QCD critical point exists and, if any, where 
it is located. 

Because of the notorious sign problem with finite quark 
chemical potential ^q, the importance sampling breaks 
down in finite-density simulations. Although theoreti- 
cal attempts are making steady progresses, (temporal) 
finite- volume effects are not easily treatable [12], and it 
is still difficult to extract any reliable conclusion even 
on a qualitative level. Then, under these circumstances, 
there are three major passages toward the QCD phase 
diagram studies (except for recent developments in the 
functional method [13]). 

1. One can discuss the critical phenomena assuming 
the QCD critical point. This is a common strategy of 
theory in general. Since the critical properties are univer- 



sal, one can make model-independent predictions. The 
virtue of this approach is the generality, but it does not 
give any clue about the concrete structure of the QCD 
phase diagram. 

2. One can utilize the effective model description 
with a reasonable choice of the model parameters [14]. 
The location of the critical point is sensitive to model 
details. Not that all model results are model dependent, 
but the nature of the phase transition at high density 
strongly depends on a part of the model setup, as we will 
elucidate later. 

3. One can make a conjecture on the phase struc- 
ture based on generic properties of QCD such as symme- 
tries [15, 16] and the degrees of freedom in a particular 
limit [2]. Because the argument lacks for concrete dy- 
namics unlike the model study, one should check individ- 
ually which scenario is favorable in reality. Nevertheless, 
such a conjecture from physics deliberation provides us 
with a useful guideline for model analysis. 

The aim of the present work is to establish a path 
from ^ to 5 in the above classification. This is a route 
rather opposite to conventional approaches. Instead of 
choosing a particular model description, we shall extract 
the essential ingredients common in most model studies 
and try to unveil the underlying physics mechanism in 
a way free from model artifacts. In particular, by look- 
ing at the saturation curve, i.e. the energy per particle 
as a function of density, we can clearly see the nature 
of the liquid-gas phase transition, which also enables us 
to understand why the vector interaction would disfavor 
the first-order phase transition. It is a straightforward 
extension to include inhomogeneity as the chiral spiral 
for simplicity, and we can then find that the phase struc- 
ture still has rich contents, which is again understandable 
from the saturation curve. 



II. FIRST-ORDER PHASE TRANSITION AT 
ZERO TEMPERATURE 

Let us start our analysis utilizing the same setup as 
Ref. [17]. We treat cold and dense quark matter in 



a quasi-particle description. This means that we as- 
sume a Fermi hquid of quark matter, which should be 
vahd for bulk thermodynamic quantities as long as T is 
small enough and the Landau damping is a minor effect. 
Strictly speaking, our strategy would work in a density 
region between two onsets; one for quark deconfinement 
and the other for color superconductivity. It is very hard 
to quantify deconfinement and a phenomenological study 
of the equation of state [18] implies that quark-hadron 
crossover may start around the baryon density p-g ~ 2pQ 
with the normal nuclear density po — 0.17 fm""^. The 
phase structure involving color superconductivity is more 
ambiguous and severely dependent on the models around 
Pb ~ 5/9o [19]. Therefore we should restrict the validity 
of our treatment within a range 2po ^ Pb ^ 5po- This 
is, however, a rather conservative estimate and should 
be loosed at higher temperature where quarks would be 
more liberated. 

In this way the thermodynamic potential from quasi- 
particles, ^matter, is expressed as a function of the effec- 
tive mass M in a form of 

anattcr[A/]/l/=- f^^p' p{p') 
JQ 

-m^N,Tj^\n{l + e-^/^), (1) 

where p{p) is the quark number density defined by p{p) ~ 
2NcN{ J -^;^[n-p{ujp — p) — riF('^p + m)] with the Fcrmi- 
Dirac distribution function, n-p{ujp) = (e"''/^ -I- 1)^^, and 
the quasi-particle energy, ujp ~ y^P^ + . It is im- 
portant to note that this ^q-dependent matter part is 
common in any quark models such as the (P)NJL and 
the (P)QM models [14]. Then, the model uncertainty is 
unavoidable in the vacuum part. 

In a quasi-particle picture of quarks the vacuum part 
could be expressed as fJo[A'/]/y = —2NcN{ (1^)3 + 
U[M] with a potential term. If we postulate it as 
U[M] = (M - m)V(4g,), then QoiM] + r!„,attcr[M] ex- 
actly amounts to the thermodynamic potential in the 
NJL model with the bare mass m [20]. To implement 
the U{1)a anomaly in the three-flavor case, we may add a 
term —gd{M—m)^ in [/[M]. From now on, we shall adopt 
a more general form of flo [M] inspired by the Ginzburg- 
Landau expansion, i.e. 

no[M]/V = a{M^ - M^f - bM - cM^ . (2) 

Although the thermodynamic potential in hand is ex- 
tremely simple, this setup sufficiently grasps the generic 
features of the phase transition in cold and dense quark 
matter. One may wonder if this polynomial form would 
miss a logarithmic singularity as discussed in Rcf. [21]. 
There are two reasons why this is not a serious problem 
to our analysis: First of all, such a logarithmic singularity 
is related to the infrared singularity of massless fcrmion 
loops. As we will see later, we are more interested in the 
massive case than the chiral limit and the effect of the 
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FIG. 1. Potential shapes from Eqs. (1) and (2) at T = 
with pq = 370 MeV. r^mattcr is model independent, while 
the vacuum potential f^o leaves uncertainty. The parameters 
in fio are chosen as Mo = 340 MeV and the curvature from 
a — 0.05 (« a value in the Unear a model) to a = 0.08 (~ a 
value in the NJL model), and b = c — 0. 

logarithmic singularity is only minor then. Second, this 
logarithmic term has no effect for the first-order phase 
transition at large pq and T = because the phase tran- 
sition typically exists around M ^ Mq (see Fig. 6), which 
is far from the singularity near M = 0. 

To enter the regime at higher temperature, one should 
consider the meson fluctuations that may give rise to 
T-dependcnt coefficients in Eq. (2). Therefore, strictly 
speaking, our analysis is valid only in the region with 
Pq ^ T. In what follows we consider only the c = 
case, for we are interested in the mechanism in favor of 
the first-order phase transition and c ^ would trivially 
stabilize the first-order transition. 

Figure 1 shows the typical behavior of the potential. 
As discussed in Ref. [17] the matter part fJmattor al- 
ways has a minimum at M = because the baryon 
density is the largest when quasi-particles are mass- 
less. Let us consider the condition for the first- 
order phase transition in the case of T = in which 
Eq. (1) simplifies as: flmattor/1^ = -iNcN{/l2Tr'^){pFp^ - 
^M^PF P + ^MH n[{p + pf)/{p - pf)]) 9{p - M) with 
Pp = ^ p?' — M^. In Ref. [17] the upper bound for the 
curvature a was estimated under a reasonable but lim- 
ited situation, pq ~ Mq. We can relax this numerically 
only to find that a first-order phase transition can re- 
main in the chiral limit (6 = 0) unless we choose unphys- 
ical parameters so that a phase transition takes place at 
Pq 3> Mq. Then fimattcr strctchcs far beyond M ~ Afo 
and the phase transition is no longer of first order. 

This simple analysis tells us that the first-order phase 
transition at T = can occur since fimatter is propor- 
tional to 9{p — M) and fl does not have to contain a 
term, while fl is sometimes assumed to take a form 
of C2M2 + c^M^ + cgAfS at T 7^ 0. Thus, the present 
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FIG. 2. The location of the potential minima at the first- 
order transition point when = (solid curve) and = 
0.12/Mo^ ~ fO^^MeV^^ (dashed curve) are chosen, respec- 
tively, with a = 0.05 and Mo = 340 MeV. 

formalism based on the quasi-particle approximation is 
more appropriate for the investigations of cold and dense 
quark matter. 

Furthermore, we must add a term oc in i^mattcr, 
which stems from the vector-channel interaction (ipjfj.ip)'^ 
that is chiral symmetric [22], i.e. 

17vec[M]/1^ = gvp' , (3) 

which can be evaluated with p numerically which is ob- 
tained as p = ^(/i^ - M^f'^e{pL - M) at T = 0. We 
should note that in the mean-field NJL model with the 
vector interaction, usually, the vector interaction would 
shift the chemical potential, which pushes the energy up 
by ^ 2(7vP^, and the condensation energy is negative, 
—g^p^, leading to ^ 'i.g^p^ — g^p^ = g^p^ in total. Here 
we simply postulate this in a form of Eq. (3). 

For a deeper insight. Fig. 2 is quite instructive. This 
figure shows the location of two degenerate minima in 
the potential (i.e. the dynamical mass) when /iq takes a 
value at the first-order phase transition. For example, in 
the chiral limit, the dynamical quark mass jumps from 
M ~ Mq to M = 0. The jump is naturally reduced at 
larger h (larger quark mass) and eventually only crossover 
remains beyond the bend of the curves in Fig. 2. One can 
notice that the curve substantially shrinks with positive 
f/v which disfavors the first-order phase transition. 

It is interesting to see that the vector interaction has 
only a minor impact for 6 = 0. This is because the min- 
imum at M = is intact as long as chiral symmetry is 
exact at & = and p and thus the vector interaction is 
still very small at M = Mq. This observation is, however, 
not completely free from the model choice. If the phase 
transition is located at pq > Mq with some other choice 
of parameters, the potential minimum around M = Mq 
is also influenced substantially by the density effect and 
thus the first-order phase transition could be diminished 
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FIG. 3. Energy per particle as a function of the density. The 
solid curves represent the homogeneous results for (1) b = 
gv = 0, (2) b = 0.08M|? and = 0, (3) b = O.OSAfo^ and g^ = 
0.12/A/o from the bottom to the top. The dashed curves with 
the label "CS" represent the chiral-spiral results for respective 
parameters. The horizontal axis is given in the unit of the 
normal nuclear density po = 0.17 fm~^. 

by the vector interaction. This part of uncertainty is not 
relevant, for we arc interested in the physical world with 
finite quark mass after all. 

Guided by Fig. 2 we shall specifically look at the fol- 
lowing three cases: (I) = gv = (first-order), (2) 
b = O.OSMq and gv = (weak first-order), and (3) 
b = 0.08M3 and g^ = 0.12/M^ (crossover). 

For later convenience we shall plot the energy per parti- 
cle e/pe atT = Oin Fig. 3, where £ = U/V+pbPb-^o/V 
is the internal energy density measured from the hadronic 
vacuum with M ~ Mq (before a finite density appears), 
and pb = p/Nc is the baryon number density. If the curve 
has a minimum as a function of pe, i.e. d{e / p^) / dp^ = 
tJ-B/pB — £/ Pb ~ ^1 pressure difference becomes zero, 
which indicates a first-order phase transition of the gen- 
eral liquid-gas type (see Ref. [23] for a review and also 
Ref. [24] for experimental studies). Therefore, whenever 
e/ Pb has a minimum as a function of ps , the T = sys- 
tem must have a first-order phase transition in the same 
way as the (symmetric) nuclear matter phase transition 
at /iB = Mn - B with Mjv 939 MeV being the nu- 
cleon mass and S ~ 16 MeV the nuclear binding energy. 
At the second-order transition, the energy curve should 
be flat at the point of inflection. This kind of analysis 
on quark matter is well known in the context of quark 
droplets [25] but less applied in the phase diagram re- 
search. What is necessary for the existence of the critical 
point (first-order phase transition) is a convex structure 
of the curve (saturation property), which is a general 
statement that does not rely on any model nor Ansatz. 

Because this point of the liquid-gas transition is so im- 
portant, let us recall here how an intermediate density 
between p = and the saturation density p ~ pQ can be 
realized in this case. If the energy per particle has a min- 
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FIG. 4. Schematic figure of tlie saturation curve of nuclear 
matter witli a minimum at po — 0.17 fm^'^ and tlie binding 
energy given by the volume term av in the Bethe-Weizsacker 
mass formula. An intermediate density p < po can be realized 
as a spatial average over bubbles with the core with p ^ po 
in the empty vacuum. Though the surface energy effect is 
not considered in the above schematic figure where a simple 
nucleon-gas picture is depicted, the actual bubble shapes in a 
nuclear liquid depend on the surface term as, etc. 



imum as schematically shown in the upper panel of Fig. 4 
it vifould be energetically preferable to form bubbles with 
the core with p ^ po rather than a homogeneous dis- 
tribution of dilute p. If we consider the surface energy, 
the density gradient (Weizsacker) term, and the charge 
neutrahty, bubbles should take optimal shapes such as 
the nuclear pasta (spaghetti, lasagna, etc) [26]. Such a 
state of matter is nothing but a mixed phase associated 
with the first-order phase transition, and importantly, 
this argument already implies the existence of an inho- 
mogcneous ground state near the liquid-gas transition. 
In other words, if a mixed phase is characterized by a 
typical wave number g, how can we strictly distinguish 
such a phase from an inhomogeneous ground state? One 
may think that in the case of quark matter the inhomo- 
geneity is turned on not in the density only but in the 
mass M unlike nuclear matter. We would stress, how- 
ever, that M also controls the density and the physics is 
just the same if seen in terms of the saturation curve as 
in Fig. 3. 

It is obvious from Fig. 3 that the vector interaction 
as in Eq. (3) disfavors the first-order phase transition. 
The minimum in e/pB is pushed up by the quadratic 
term cc pg and eventually the first-order phase transition 
disappears when the minimum is lost, as demonstrated 
by three solid curves in Fig. 3. In the chiral limit b ~ 
the branch of M = is separate, so that the first- 
order phase transition survives regardless of the vector 
interaction, which may change with different parameters 
as we already pointed out. With finite b, however, two 
branches with small and large M are smoothly connected 
and the minimum diminishes for large b and in accord 
to Fig. 2. 



III. CHIRAL SPIRALS 

One may find the usefuhiess of the saturation curve 
for analyses with a wider range of model space. From 
now on we shall consider the possibility to form inhomo- 
geneous chiral condensates. We here utilize the simplest 
Ansatz to introduce it, namely, the one-dimensional chi- 
ral spiral; ('i/'V') = X'^os{2qz) and (tpj^Tsip) = xsin(2(72;) 
(see Ref. [27] for reviews). This ground state of the chi- 
ral spiral can be equivalently described by a chiral ro- 
tation ip = e^'^'^'^^'^^V'' with a homogeneous condensate 
^ = {ip'ip') in the chiral limit. Then, the quasi-particle 
dispersion relation in the ■(/''-basis is expressed as [27, 28] 



(4) 



where ± in front of q corresponds to the flavor and the 
chirality that also depends on the sign of Pz ■ 

This type of inhomogeneity pattern has been consid- 
ered repeatedly in various contexts such as the pion con- 
densation in nuclear matter [28], large- A'c QCD [29], the 
Overhauser instability [30], the quarkyonic spiral with 
confining force [31], and so on. The dispersion rela- 
tion (4) should be plugged into r^matter/^ in Eq. (1). 
Unhke the normal dispersion relation, we see that a large 
part of the mass effect can be absorbed by g ~ M, with 
which p is no longer suppressed even at large M. This is 
the reason why a first-order phase transition can occur 
from the homogeneous hadronic phase to the chiral spiral 
where M is substantially large. Also, we should point out 
that the Ginzburg-Landau analysis in Ref. [32] to con- 
clude that the chiral spiral is less favored might be inad- 
equate; the largest energy gain in fl matter /V comes from 
the region with large M where the Ginzburg-Landau ex- 
pansion should not work. 

The physical mechanism to lower the total energy is 
the Overhauser effect as argued in Ref. [30] . In the ordi- 
nary Overhauser instability the momenta of the spin-up 
component are shifted up by pf and those of the spin- 
down component are shifted down by pp, so that a gap 
opens where two energy dispersion relations cross. In 
(l-l-l)-dimensional NJL model the situation is completely 
analogous [27]; a choice of g = 2/iq eliminates the p-q de- 
pendence and the energy gain originates from the fact 
that p is completely insensitive to M and thus p is never 
suppressed by M in contrast to the homogeneous solu- 
tion. In (3+l)-dimensional case, on the other hand, not 
only Pz but also p± share the Fermi momentum, and so 
the optimal q is not 2pq but rather q ^ M which will be 
confirmed by numerical calculations later. 

Thus, Oniattor always tends to favor the chiral spiral 
with q ^ M, while it is f^o that would hinder the growth 
of q. In the leading order the vacuum part has an expan- 
sion in terms of q as 

no[M, q]/V = no[M, q = 0]/V + [aM^ + pb)q^ , (5) 

where the first term with a > is a "kinetic" term 
against spatial modulation. This term should be van- 
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FIG. 5. Typical phase diagrams with the chiral spiral. The 
solid curve in the lower-/iq side represents the homogeneous 
chiral phase transition of first order with fo = = with 
which the chiral spiral region surrounded by the first-order 
phase boundaries is attached. For b = O.OSMq and (?v = 
12 1 Ml the homo eeneous first-order transition and thus the 
QCD critical point no longer appear, but the inhomogeneous 
region is enlarged as shown in the higher-^q side with a first- 
order boundary (solid curve) terminating at P followed by a 
second-order boundary (dashed curve) . 



ishing at cither A/ = or g = 0, so the expansion 
should start with M^g^. One can estimate a using a 
chiral model, but one should be careful not to pick an 
unphysical term ~ lS?q^ up from gauge- variant regular- 
ization. The latter term oc /3 comes from a phase of the 
current mass term associated with the basis change from 
"0 to V''- Quantitative details may depend on a and /3, 
but qualitative features as we discuss below do not rely 
on any specific choice of them. 

Figure 5 shows typical behavior of the phase bound- 
aries on the /iq-T plane with zero and non-zero h and g^. 
For demonstration we chose a = 0.25 and /3 = 0.25/Mo. 
Then in the lower-^q side of Fig. 5 we see that there 
is an island structure of the chiral spiral surrounded by 
the first-order boundaries. The solid curve extending 
to smaller /.tq represents a first-order phase transition 
associated with the homogeneous condensate only. It 
should be mentioned that the first-order phase transi- 
tion at & = in the high-T and small-/iq region, which 
is not of our present interest, might have been artifi- 
cially strengthened due to the lack of the logarithmic 
singularity in Eq. (2). The first-order boundary of inho- 
mogcncity at smaller yUq stays very close to this curve. 
This is because the effective potential becomes very shal- 
low near the first-order phase transition in the homo- 
geneous case as clearly recognized in the total potential 
presented in Fig. 1. The secondary first-order boundary 
at larger /iq is much weaker because M and thus q are 
small there. (Note that, in the chiral limit, q may not 
decrease but only increase in a narrow region of ^q as 
shown in Ref. [30]. This tendency near the first-order 



phase transition is partially seen also in the massive case 
in Fig. 7.) The corresponding saturation curve of e/pB 
is shown by a long-dashed curve with the label "CS" in 
Fig. 3, from which a minimum at lower energy is ap- 
parent. We note that the inliomogeneity island in the 
vicinity of the first-order phase transition is consistent 
with our intuitive discussions of the mix phase formation 
below Fig. 4. 

With the vector interaction included, the so-called 
QCD critical point is easily washed out [17, 32, 33]. Inter- 
estingly, however, as shown in the higher-^q side of Fig. 5 
and especially at P in this figure, there is a chance that 
the critical point (strictly speaking, tri-critical point) is 
revived driven by the inhomogeneous condensate. The 
question is then how robust this observation is. In fact it 
has been reported that the soliton solution [34] is more 
stable than the chiral spiral and also it exhibits a second- 
order phase transition rather than a first-order one [32]. 

Let us then consider when the second-order phase 
transition is possible in view of the saturation curve in 
Fig. 3. To have a second-order phase transition from the 
hadronic phase (with homogeneous M ^ Mq) to a gen- 
eral inhomogeneous state, there must be an energy curve 
that is tangent to the hadronic branch (solid curves from 
Pb = 0) and going below it. The curves do not have to be 
flat because there is a small energy difference before and 
after a finite density appears, which is further enhanced 
by 1 1 p\ in the slope of the saturation curve. To avoid a 
first-order transition, moreover, the energy curve should 
be monotonically increasing with increasing pB ■ 

Such a situation is not allowed, for example, in the far 
bottom curves (at b = gv ~ 0) in Fig. 3. In this case 
with the saturation energy lower than that at pB = 0, 
we can conclude that only a first-order phase transition 
is possible however complicated and optimized modula- 
tions we introduce. The situation is different with finite b 
and/or g^- It is clear on a qualitative level that a larger 
g^ would ease better inhomogeneous states to develop, 
for the dashed chiral-spiral curve could be then easily 
extended down to pB = monotonically. This means 
that the phase transition between the homogeneous and 
inhomogeneous states could be of second order. There- 
fore, unfortunately, the existence of the critical point P is 
again not a robust conclusion especially with the vector 
interaction. 

From a plain physical interpretation, it would be the 
most natural to have continuous phase transitions that 
border the inhomogeneous island. Such an intuition is 
based on the picture of the liquid-gas phase transition. 
In fact, if the boundary is a first-order phase transition, 
there will appear a density regime that can be described 
only as a mixed state. It is the role of the chiral conden- 
sate in quark matter that makes a difference from the sit- 
uation in nuclear mater. The density modulation inher- 
ent in a mixed state can be mimicked by the modulation 
in the chiral condensate, which would lead to an inhomo- 
geneous ground state of quark matter with lower energy. 
This is exactly what happens with the soliton solution 
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FIG. 6. Behavior of the constituent mass M as a function of 
and T in the unit of Mo in the case with b = 0.08A/(f and 
Qv = 0.12/Mo . At the phase boundary into the chiral spiral, 
M drops but still remains finite, and goes smaller continuously 
with increasing /iq and/or T. 
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FIG. 7. Behavior of the wave-number g as a function of /iq 
and T in the unit of A'Iq in the case with b — O.OSMq and 
gv = 0.12/Mo. This clearly shows the structure of the chiral- 
spiral island, surrounded by a steep "cliff" at smaller /iq and 
a gentle "beach" at larger /iq. 

in Refs. [32, 34]. Indeed, at the onset of solitonic inho- 
mogeneity, localized domain- walls start appearing, which 
approaches sinusoidal patterns at larger fj.q- The density 
profile has peaks arising from the kinks and this situation 
is reminiscent of a mixed state picture as schematically 
depicted in the bottom of Fig. 4. It would interesting 
to figure out the saturation curve corresponding to the 
solitonic solution. This is beyond our current scope, but 
it presumably goes below the chiral-spiral curves and is 
smoothly merged with the hadronic branch at smaller pB- 
Finally let us take a closer look at the solution with 
b = 0.08M^ and = 0.12/M^. Figures 6 and 7 show 
the behavior of the constituent mass M and the wave- 



number q, respectively, as functions of /iq and T. At a 
glance one may notice that q ^ M holds in the chiral- 
spiral region as we discussed. The structure of the chiral- 
spiral island is quite characteristic. In view of Fig. 7, one 
might say that the island is surrounded by a steep "cliff" 
at smaller /iq and a gentle "beach" at larger /iq [39]. 

Such a structure of the island should be quite robust 
because the energy gain is mainly attributed to g '~ M . 
Hence, the cliff stands with a large energy gain at smaller 
/iq where M is still large, and the inhomogeneous state 
gradually becomes indistinguishable from the homoge- 
neous state as M gets smaller at larger /iq. 



IV. SUMMARY 

We have developed a picture of the first-order phase 
transition of quark matter based on the saturation curve 
and the liquid-gas phase transition. From this picture 
we discuss the relation between the order of the phase 
transition and the behavior of the saturation curve. We 
demonstrated this using a simple Ansatz of the chiral spi- 
ral, but the argument itself is not limited to such a special 
choice. As a matter of fact, because the chiral spiral can 
be mapped to the conventional pion condensation [28] 
that is killed by the spin-isospin interaction, it may be 
likely that the chiral spiral should be suppressed by the 
axial- vector interactions ~ (V'757p'0)^ or ~ ("0757/^'^^')^, 
and eventually superseded by others such as the soliton- 
like modulation and more generally multiple- wave super- 
positions. Even in this case the saturation curve would 
provide us with valuable information on the nature of the 
phase transition. 

We can think of several directions as future extensions. 
It may be interesting to seek for some connections be- 
tween our saturation considerations and the Ginzburg- 
Landau analyses as in Ref. [36]. Also, the interplay be- 
tween the chiral spiral and the external magnetic field 
would deserve further investigations [37]. We are actu- 
ally working in this direction to clarify the phase struc- 
ture with three axes, /iq, T, and B, including the spatially 
inhomogeneous state [38]. 
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